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I present a selection of recent lattice data by major collaborations for the pseudo-Goldstone boson masses in 
full {Nf = 2) QCD, where the valence quarks are chosen exactly degenerate with the sea quarks. At least the more 
chiral points should be consistent with Chiral Perturbation Theory for the latter to be useful in extrapolating to 
physical masses. Perspectives to reliably determine NLO Gasser-Leutwyler coefficients are discussed. 



1. INTRODUCTION 

In lattice QCD physical observables like masses 
and decay constants are extracted from the fall- 
off pattern of correlation functions; for instance 
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with Q{x) = u{x)Td{x) and F G {7a 17475} allows 
to determine and F^. Like in the continuum, 
the fermions are integrated out and (|l|) reads 
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where Tr goes over colour and spinor indices. In 
this form, the quark masses msca, fn^ai should be 
equal, since they stem from the very same term 
Sf=Jdx q{D+m)q in the original form, eqn. 

However, state-of-the-art algorithms for the de- 
terminant det(D-l-m) and the Green's function 
S{x, 0) = ((-D-l-TOvai)^^)£cO get parametrically slow 
as the quark mass is taken light. Since this effect 
is much more pronounced for the determinant, 
most groups are restricted to a regime where the 
sea-quark mass for the two light ffavours (the one 
in the determinant) is about half of the physical 
strange mass, while the valence-qnark mass for 
the u/d-quarks (the one in the propagator) may 
be pushed down to '^1/4 of the physical . 

With these data at hand, one has to extrapolate 
both in TOgoa and TOvai down to the physical u- and 
d-masses (which we take degenerate); see Figl. 



Since the extrapolation range is large, it is cru- 
cial to know the functional form against which the 
data may be matched. Partially quenched chiral 
perturbation theory (PQ-XPT) provides such a 
framework 0|, but its limited range of applicabi- 
lity (around the twofold chiral limit) might cause 
a problem; it is a priori not clear whether the bulk 
of the existing data lies in this range. 

2. ELEMENTS OF XPT 

The mandatory test is most conveniently done 
along the diagonal line in Fig. 1, which is the do- 
main of (ordinary) XPT PP]. This is the effective 
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Figure 1. 
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Schematic view of partially quenched data 
taking, with crosses indicating measurements in the 
(jTisea, myai) plane for degenerate u/d-quarks. While 
the real-world pion (*) is on the unitary line (diago- 
nal) , substantial savings in CPU-time stem from tak- 
ing the detour through the sector with mvai < rnsca,- 
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theory of QCD at energies small compared to the 
intrinsic scale, 47rf,r, which faithfully reproduces 
the structure in the chiral limit. Its Green's func- 
tions are dominated by pseudo-Goldstone contri- 
butions and computed in a systematic expansion 
in the external momenta and the quark mass, the 
latter counted as m~p^ (|]. 

The leading part of the Lagrangian, L^^\ con- 
tains two low-energy constants {B, F) which are 
finite dimensionful quantities^ and the quark 
mass matrix M = diag(mj;, to„ [, mj]). The 0{p^) 
suppressed part, L^^^ contains about a dozen 
so-called Gasser-Leutwyler (GL) coefficients, de- 
noted li for 2 flavours |^ and Li if the strange 
quark is included . In the dimensionally regula- 
ted theory, the li or Li consist of a divergent part 
(cxe~^, where e = 4 — d) and a finite part (cx e") 
which is scale-dependent. After some reshuffling 
of finite contributions in MS style, the latter part 
(denoted by or L\{ii)) encodes for the QCD 
short distance contributions. In an arbitrary 
Green's function calculated at order O(p^), the 
divergence from the loop compensates exactly the 
divergent part of the GL-coefficient (see Fig. 2), 
and an important check of the finite part is that 
the scale-dependent (finite) integral J(/z) and the 
renormalized "constants" (or L\{^)) add up 

to form an exactly scale-independent expression. 

What we are interested in is the prediction for 
the squared pseudo-Goldstone boson mass in 2- 
flavour QCD versus the sum of the degenerate 
valence quark masses, i.e. Af^ versus 2m. For the 
latter, we adopt (MS,2GeV) conventions, which 
implies that the leading-order (LO) constant B is 
also quoted in MS scheme at ^ ~ 2 GeV, since the 
product mB is scheme- and scale-independent (as 
is F). An important simplification stems from 



is the pion decay constant in the chiral limit ^ 87 MeV, 
while _B = — E/F^ with E the chiral (one-flavour) conden- 
sate (in the 2- or 3-flavour theory) in the chiral limit. 
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Figure 2. Contributions to the pion self-energy at 
NLO: 1-loop graph with a vertex from L*-^^ (little dot) 
and a counterterm from i''*' (fat box). The divergent 
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compensate each other, and in the finite 



the fact that the GL-coefhcients l\{^JL) may be 
traded for the low-energy scales A^ , and in terms 
of these also the prediction at next-to-leading or- 
der (NLO) @ in XPT takes a form in which the 
scale-independence of the physical result is ex- 
plicit (see jllP for details and references): 
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3. AIMING AT NLO GL COEFFICIENTS 

Equations (|, |) may be used in two ways: 
Either phenomenological values for A3 are 
plugged in and the resulting curves are compared 
to lattice data, or the functional forms are fitted 
to the data, which amounts to a lattice determi- 
nation of the QCD low-energy constants. 

The first option is shown in Fig. 3, where quark 
masses determined via the axial (AWI) or vector 
(VWI) Ward-Takahashi identity |§ are renor- 
malized by 1-loop perturbation theory |^,||. The 
heaviest pion is GeV, and one sees that it is 
still compatible with the LO graph, but not with 
the NLO or NNLO version. The reason behind 
is that XPT yields an asymptotic series: the LO 
result for may be good up to ^1 GeV, while 
the NLO/NNLO forms are more precise at small 
quark masses, but break down earlier, e.g. the 
NLO functional form might be good up to (say) 
600 MeV, to NNLO form up to (say) 400 MeV. 



CP-PACS VWI, pert. ren.. 
CP-PACS AWI, pert. ren.. 
UKQCD VWI, pert. ren.. 
UKQCD AWI, pert. ren.. 
LO-XPT (B from phenomenology) 
NLQ-XPT (B F Ao from phenom. 
NNLO-XPT (A„ from pfien, k„=oS 
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parts (oc the /i-dependence cancels exactly. 



Figure 3. LO/NLO/NNLO chiral predictions for 
A/^ (with phenomenological values for the low-energy 
constants) compared to CP-PACS and UKQCD data, 
using the VWI and AWI definition of the quark mass. 
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CP-PACS VWI, pert. ren. 
CP-PACS AWI,_pert. ren. 

B F A, fitted 
B A3 fined, F fixed 
UKQCD VWI, pert. ren. 
UKQCD AWI.pert. ren. 

B F A, fitted 
B A, fitted, F fixed 
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Figure 4. Fits of the NLO functional form (g) to the 
CP-PACS |] and UKQCD [|| data, after dividing by 
2m. Either F is kept fixed with the 3 most chiral 
points used or B, F, A3 are fitted, using all 4 points. 

The second option is likely premature: If the 
estimate that the NLO functional form is good 
up to (M^rro)^ ~2 is correct ^ (with rg ~0.5 fm), 
then Fig. 1 depicts the current situation, and 
there is no room for using eqn. (|^) to fit the data, 
since it has 3 parameters and in cither set at best 
1 data-point is in the permissible range. Still, one 
might simply go ahead and see how inconsistent 
results get, if one ignores that the bulk of the data 
is likely not in a regime where (PQ-)XPT applies. 

For practical reasons, a factor 2m is taken 
out and parameters are made dimensionless with 
powers of rg, i.e. the actual NLO form used is 



2mrQ 



^Bro 



log- 



(5) 



with the pertinent fits shown in Fig. 4 and the 
resulting parameters summarized in Tab. 1 

Clearly, results should be taken with care, but 
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Table 1 

Coefficients in the fits of the functional form 
the degenerate CP-PACS § and UKQCD [|] data 
with average VWI and AWI masses j^. Constrained 
values in brackets, ranges such that ^^(d.o.f.)^ 1. 



it is reassuring to see that the agreement among 
the fitted values of the NLO GL constant A3 is not 
substantially worse than for the LO constant B. 
The numerical values in the last column of Tab. 1 
correspond to A3 ~ 1.4 GeV. This is consistent 
with the early estimate A3 ~ 0.6 ± I I GeV §, 
and hence supports the hypothesis that at least 
in two-flavour QCD the chiral condensate is large 
(see Q for references to this topic). However, 
the difference between f ^ from the unconstrained 
fit and the physical value (cf. the central column 
in Tab. 1) reminds one that to date the bulk of 
the data is likely outside the regime of validity 
of (PQ-)XPT and accurate simulation results at 
smaller quark masses are needed (cf. Q) before 
one may draw definite conclusions. 

In summary, vs. 2m illustrates that lat- 
tice QCD is, in principle, an ideal tool to deter- 
mine certain NLO GL coefficients. To date, most 
studies with Nf =2 Wilson quarks are likely in 
a mass-range where there is only a "point-like" 
overlap with the regime where (PQ-)XPT holds. 
Once moderately lower masses are covered and 
with further improvements in renormalizing the 
bare data and controlling the lattice artefacts (see 
1^ for a discussion) the lattice will yield a compe- 
titive determination of those NLO GL coefficients 
which encode how low-energy QCD Green's func- 
tions depend on the quark mass. 
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